Signatures of the superfiuid to Mott insulator transition in cold bosonic atoms in a 
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We study the Bose-Hubbard model using the finite size density matrix renormalization group 
method. We obtain for the first time a complete phase diagram for a system in the presence of 
a harmonic trap and compare it with that of the homogeneous system. To realize the transition 
from the superfiuid to the Mott insulator phase we investigate different experimental signatures of 
these phases in quantities such as momentum distribution, visibility, condensate fraction and the 
total number of bosons at a particular density. The relationships between the various experimental 
signatures and the phase diagram are highlighted. 
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I. INTRODUCTION 



In recent years many theoretical and experimental in- 
vestigations have been carried out in the field of ultra- 
cold atoms in optical lattices [J [3, Q. The interest 
in bosonic systems began with the seminal paper by 
Fisher et al where a phase transition from a su- 
perfiuid (SF) to a Mott insulator (MI) in a lattice of 
bosons was predicted when the on-site Coulomb repul- 
sion between the atoms dominates the nearest neighbor 
hopping amplitude. Since then, a variety of theoretical 
approaches || 0, 0, S, ©, M, M, [13, El HI El have 
been used to study the Bose-Hubbard Model (BHM) 0. 
There is good agreement between the phase diagrams ob- 
tained from the different techniques. While the BHM was 
originally developed in the context of 4 He Q , its poten- 
tial to describe ultra cold bosons trapped in an optical 
lattice was soon realized by Jaksch et al [lr3 | . This paper 
has had a great impact on the condensed matter commu- 
nity because high-precision experiments on cold atoms in 
traps can now be used as a powerful and reliable tool to 
study a variety of quantum phase transitions [l|, 0, 9- 
The experimental realization of the quantum phase tran- 
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sition from the superfiuid to the Mott insulator in three 
dimensio ns Jl7l| , two dimensions [Tc| , as well as in one di- 
mension |19( | soon followed. The bosons in an optical lat- 
tice are well described by the Bose-Hubbard model mod- 
ified to include a trap potential 16], which is normally 
harmonic. In the presence of a trap, the density profile 
exhibits a rich structure as the SF and the MI phases co- 
exist. A variety of numerical methods have been applied 
to understand the model [H HE HI HI HI HE Hf . The 
most important aspect that has emerged is the lack of a 
global phase in these systems. As an analogy, in three 
dimensions the superfiuid and the Mott insulator phases 
coexist as shells in an onion. The unprecedented control 
over the system parameters by tuning the laser intensity 
has paved the path for the experimental realization of 
these predictions (271,128}. 

In this paper, we re-visit the one dimensional Bose- 
Hubbard model using the density matrix renormalization 
group method. Our main motivation is to obtain the 
phase diagram, given the experimental realization of the 
shell structure and their signatures. The Bose-Hubbard 
model, describing bosons in an optical lattice, is given by 
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where t is the hopping amplitude, U is the on-site re- 
pulsion between the atoms and Vt is the depth of the 
external trapping potential, a\(ai) are the bosonic cre- 
ation (annihilation) operator, n., = aja, is the number 
operator and r, is the position of the i th lattice site from 



2 



the trap center. For simplicity, the energy is scaled in 
units of the hopping amplitude which is set to one. 

In experiments, the optical lattice potential (Vol), 
formed by the superposition of three counter propagating 
laser beams [2J, |29[, can be written as: 

Vol = V Q {sin 2 {k L x) + sin 2 (fc L y) + sin 2 (fc L z)) (2) 

where Vq is the lattice depth measured in units of the 
single photon recoil energy Er, fe& = 2it/\l is the wave 
vector and Al is the laser wavelength. By increasing 
the intensity of the beam in two directions, the hopping 
is restricted to one dimension and this results in a one 
dimensional lattice. The on-site interaction U and the 
hopping amplitude t are related to Vq and Er as follows: 



j (xexp{2V /E R ) 1 / 2 . 



(3) 



Therefore, by varying the lattice depth the ratio of U/t 
can be tuned. 

Earlier studies of the homogeneous system show that 
when the on-site interaction strength U is small com- 
pared to the hopping amplitude t, the system remains 
in the SF phase characterized by long range coherence. 
When U increases and becomes much larger than t, a 
transition from the SF to the MI phase occurs at some 
critical value of U = U c ~ 3.4 (in units of t) @, S, GU- 
This transition belongs to the Kosterlitz-Thouless uni- 
versality class [10, HH . The SF-MI transiti on g ets more 
interesting in the presence of a trap [H, l33l. 1341 [H[. 
In this case, the entire system remains in the SF phase 
for small U values, but as U increases, a MI phase de- 
velops around a central SF phase, followed by a SF 
shoulder. Further increase in U ultimately results in 
a MI phase thoughout the lattice with the exception 
of the trap edges, which we refer to as a central MI 
phase. This alternating occurrence of the SF and the 
MI phases can be observed in the number density profile. 
The system exhibits a plateau at integer densities (MI 
phase) surrounded by a region of non-i nteg er densities 
(SF phase) [U M, El H Hi HE US M, [13 Hpland 
this has been recently observed in experiments \2ii HH in 
3D optical lattices. 

In this paper, we first obtain the phase diagram of the 
model, given by eqn. ([1]), in the homogeneous (Vt = 0) 
limit. We then extend our analysis to the inhomogeneous 
case and compare it with the homogeneous case. Finally 
we obtain experimentally measurable quantities like con- 
densate fraction, visibility and density profile and deduce 
the inhomogeneous phase diagram. We also discuss how 
the shell structure in the optical lattice can be observed 
from condensate fraction and visibility. 

This paper is structured as follows: Section (|TT| con- 
tains a brief discussion on the Finite Size Density Ma- 
trix Renormalization Group technique (FS-DMRG). Sec- 
tion (jllip discusses the phase diagram obtained from the 
number density profile, while Section (|IVp analyzes the 
experimental signature for the phase transition and we 
summarize our results in Section (fVT). 



II. METHOD OF CALCULATION 

We have employed the finite size density matrix renor- 
malization group method (FS-DMRG) [H, Hz[ with 
open-boundary condition to determine the ground state 
energy and the wave function of the system. This method 
is one of the most powerful technique in one dimension 
and has been widely used to study the Bose-Hubbard 
model @, i, 0, HI (H- We have considered a soft- 
core case by retaining four bosonic states per site and 
the weight of the states neglected in the density matrix 
formed from the left or right blocks is less than 10~ 6 . For 
better convergence of the ground state energies of various 
phases, we have performed a finite size sweeping proce- 
dure [|, [H, twice in each iteration of the FS-DMRG 
method. 

Our FS-DMRG method consists of two steps: (i) 
a DMRG iteration where the length of the system L 
is increased to L + 2 and (ii) finite size sweeping to 
achieve better convergence of the ground state energy 
El{N) [1, H||. We consider a system with an initial 
length L = 4 and the number of bosons N — 4. The FS- 
DMRG procedure is then employed, keeping the density 
of bosons fixed at p = 1 , until we have a desired number 
of bosons in the system, say N = 30. Then onwards, in 
the FS-DMRG iteration, only the length of the system 
is increased keeping N fixed, until the system grows to a 
desired size, for example L = 100. At this point we have 
a system of length L = 100 with N = 30 bosons. Now 
keeping the length fixed, we increase N in steps of 1 at 
the end of each FS-DMRG sweep. In our example, N is 
varied from 30 to 150. Thus at the end of the FS-DMRG 
calculation, for a given set of parameter values, we have 
ground state energies El(N) and wave functions \iPln) 
for a system of length L — 100 with N varying from 30 
to 150. From El(N), we obtain the chemical potential /i 
of the system 
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The on-site local number density (rii) and the local com- 
pressibility Ki are defined as 



(Ui) = (lpLN\ni\lpLN) 
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(5) 



(6) 



For the homogeneous case, the local density (m) — 
N / L and hence the compressibility are uniform through- 
out the system except at the edges due to boundary ef- 
fects. The chemical potential corresponds to that of the 
whole system. However, for finite Vt the lattice is inho- 
mogeneous and as a result the local density (n^) attains 
its maximum value at the center of the trap where the 
potential is minimum and decreases as we move away 
from the center, eventually going to zero. 
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In most of our calculations, we have taken Vt = 0.004 
and 0.008, L = 200, with N ranging between 30 and 150 
and U between 2 and 20. 



III. RESULTS AND DISCUSSIONS 



A. Homogeneous case 




FIG. 1: Variation of the compressibility n (lower panels) and 
the density (n) (upper panels) as a function of p for the ho- 
mogeneous case (Vt = 0). The plateau region at (n) = 1 and 
k = signals the incompressible (p — 1) MI phase. 
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FIG. 2: Phase diagram for the homogeneous Bose-Hubbard 
model. The MI phase has density p = 1. 



We begin our discussions with the homogeneous case. 
The density of the system (n) as a function of the chem- 
ical potential p for three values of U = 3,4 and 7 are 
shown in the top panels of Fig. ([I}. The formation of 
a plateau in the (n) versus p plots, for U — 4 and 7 at 



(n) = 1, in contrast to U = 3, signals the onset of the 
p = 1 MI phase, where p denotes the density per site. 
The lower panels in Fig. ([1]) show the compressibility n 
(calculated using Eq. ©) as a function of p. It is clear 
that the MI phase is incompressible, i.e., k = 0. The 
cusp in k as (n) approaches 1, shown in Fig. {!]), is due 
to quantum criticality. The phase diagram for the homo- 
geneous system is obtained by picking out the values of p 
at the knees where (n) = 1 and k = and plotting them 
in the p — U plane. This is shown in the fig. ^ and is in 
agreement with earlier results p| [^, [T3, Il3l |. The cusp in 
the compressibility reflects the Kosterlitz-Thouless type 
behavior of the SF-MI transition. 



B. Inhomogeneous case 

We now analyze the case when there is a finite trapping 
potential V t in the Bose-Hubbard Hamiltonian. Taking 
the depth of the trap Vt = 0.004, number of bosons 
N = 100 and length L = 200, we obtain the local density 
profile (rij) as a function of the distance from the center 
of the trap r-j, as shown in Fig. ©. In contrast to the 
homogeneous case, (rii) is not uniform when the trap is 
finite. It decreases monotonically as we move from the 
center of the trap towards the edges. For larger values 
of U (U > 6) the density profile develops a well defined 
plateau at (m) = 1 and the length of the plateau grows 
as U is increased further. A simple way to understand 
the behavior of (rii) is through the Local Density Ap- 
proximation (LDA) [23|, HH where the local density at 
site i in the trapped case is given by the density of a 
homogeneous system with a chemical potential 



Mi 



(7) 



Here p refers to the chemical potential at the center of 
the lattice where the trap potential is zero. Using Eq. (|7|) 
we can re-scale the x-axis of Fig. ([3]) , so that we get the 
local density as a function of the local chemical potential, 
as shown in Fig. (j4|) for U = 3 and 9, where the homo- 
geneous case is also documented for comparison. The 
local density as a function of the local chemical potential 
exhibits trends similar to the homogeneous system, that 
is the value of (m) increases smoothly as pi increases 
for lower values of U , where the entire system is in the 
compressible supcrfluid phase. However, as U increases 
further, a clear plateau emerges at (rii) = 1 indicating 
the onset of the p = 1 MI phase. 

The local compressibility can be obtained from the lo- 
cal density using eqn. ©. Fig. §5§ shows the local density 
and the local compressibility in the upper and the lower 
panels, respectively. The plateau region at (m) = 1 (up- 
per panels) has the corresponding local compressibility 
equal to zero. This verifies that the plateau present in 
the number density profile represents the MI phase. Fur- 
thermore, this confirms the coexistence of the SF and 
the MI phases in the presence of a trap. Though there 
is an overall agreement between the homogeneous and 
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the inhomogeneous cases, we note that there are slight 
discrepancies. The sharp SF-MI transition observed in 
the homogeneous system is smoothened out in the pres- 
ence of a trap. The cusp like behavior observed in k 
for the homogeneous case is also lost. The agreement 
between the homogeneous and the inhomogeneous cases 
prompted us to obtain the phase diagram for the inhomo- 
geneous system, by analogously picking out the values of 
/i at the knees where (n-) = 1, and plotting them in the 
p — U plane. The resultant phase diagram is compared 
with the homogeneous result in Fig. (O . It is interesting 
to note that the MI lobe for the inhomogeneous system 
lies within that of the homogeneous system. However, 
the MI lobe retains its cusp like shape. Thus the SF-MI 
transition for density p = 1 and Vt > does have the 
Kosterlitz-Thouless universality class behavior. 




FIG. 3: Density profile as a function of n for Vt = 0.004, 
L = 200 and TV = 100 for a range of U. Note that as the on- 
site repulsion increases, a MI phase forms around the central 
SF phase, finally leading to a central MI phase for higher 
values of U. 



IV. EXPERIMENTAL SIGNATURES 
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FIG. 4: Density as a function of chemical potential for the 
homogeneous and inhomogeneous cases. The solid line repre- 
sents the homogeneous system, the dashed lines and dashed- 
dot lines, the inhomogeneous case (Vt = 0.004 and Vt = 0.008 
respectively). The homogeneous case has sharp transitions 
compared to the finite trap case. 
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FIG. 5: Local density and local compressibility as a function 
of the local chemical potential for Vt = 0.004. 



In the earlier section we have established that the 
ground state phases for the Bose-Hubbard Hamiltonian 
given by eqn. (JTJ) are either the superfluid or the Mott in- 
sulator depending on the ratio oiU/t and the local chem- 
ical potential /x, . Since ^ is uniform for the homogeneous 
system, the ground state is global in nature. However, 
for inhomogeneous systems, ^ being non uniform, both 
the SF and the MI phases coexist, as already discussed. 
It would be worthwhile to explore the signatures of these 
coexisting phases in experimentally determined quanti- 
ties. 

It is now possible in experiments to record the spa- 
tial distribution of the lattice with different filling fac- 
tors using spatially selective microwave transitions and 



spin-changing collisions as shown by Foiling et al. [271 ]. 
Similar experiments in one-dimensional optical lattices 
can yield density profiles from which the phase diagram 
can be obtained, as discussed in the previous section. An- 
other way to obtain direct information about the Mott 
plateaus (shells in 3D) is through the atomic clock shift 
experiment [28| . By using density dependent transition 
frequency shifts, sites with different occupation can be 
spectroscopically distinguished, thus giving us the infor- 
mation about the number of sites with a given density 
p, defined as N(p). In Fig. ([7]) we plot N(p) versus p 
for several values of U for a system with V t = 0.004 and 
N = 100. The peak in N(p) at p = 1 for U > 6.0 is a di- 
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FIG. 6: Homogeneous and inhomogeneous phase diagram 



rect signature of a well developed Mott insulator plateau 
in the inhomogeneous system. The size of this peak in- 
creases with U and is consistent with the increase in the 
length of the MI plateau as shown in the phase diagram 
(see Fig. ©). From N(p) we obtain the following quan- 
tities: Nmi/N, the fraction of number of bosons in the 
Mott plateau and N$f/N, the fraction of bosons in the 
SF region. Here Nmi and N$p are the number of bosons 
in the MI and the SF phases respectively. Fig. (jHJ) shows 
both Ni/n/N and N$f/N for several values of U keeping 
V t = 0.004. For U < 6.0, we observe that N SF /N is close 
to one, while Nmi/N is close to zero. This is because the 
entire system is in the SF phase. However, for U > 6.0, 
the increase in Nmi/N, signals the formation of a MI 
plateau in the system. The critical value of U, marking 
the transition to the MI phase (p = 1) can be read-off 
from Fig. ((SJ) and is given by Uc ~ 6. The small plateaus 
seen in Nmi/N and N$f/N are indicative of the detailed 
distribution of the bosons as U increases in the presence 
of a trap. 

In other experiments, the cold atom gas trapped in an 
optical lattice is allowed to expand and the interference 
pattern in the density of the expanding gas is recorded. 
The density distribution is mirrored in the momentum 
distribution defined as, 

1 L 

n ^ = T K°i) ex P(*9( fc " 0) ( 8 ) 
k,l=l 

where k, I are the lattice sites. Fig. shows the momen- 
tum distribution for different U values. The superfluid 
phase that has long-range coherence exhibits sharp in- 
terference peaks, while the Mott insulator phase, where 
the local density is pinned to integer values per site, 
breaks this coherence and hence no sharp peaks are ob- 
served (l7j . Although the presence of the interference 
peaks in the density distribution (or analogously momen- 
tum distribution) was originally used to signal the forma- 
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FIG. 7: N(p) versus p for Vt = 0.004. A small offset is added 
to the Y-axis for clarity. The peak in N(p) at p = 1 signals 
the MI plateau in density profile (Fig. [3j and thus confirms 
the coexistence of SF and MI phases. 




U/t 



FIG. 8: Nmi/N and N S f/N as function of U. The critical on- 
site interaction Uc ~ 6.0 for the SF-MI transition for p = 1 
can be easily read-off. 



tion of a SF phase, recently it has been established that 
the visibility of the interference fringes [1^, [H, GjjJ S3| 
provides a clear signature of the transition. The fringe 
visibility is defined as 

_ AT 4- N ■ 

1 'max \ 1 v mm 

where N max and A^ m i n are the maximum and the min- 
imum of the momentum distribution measured at q = 
±27r and q = ±7r respectively, in one dimension. The 
condensate fraction, that is the number of bosons in the 
condensate with respect to the total number of bosons, 
is defined as the largest eigenvalue of the matrix (ajaj) 
divided by the total number of bosons [42j]. The fringe 
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visibility and the corresponding condensate fraction for 
V t = 0.004, N = 100 and V t = 0.008, N = 50 are given 
in Figs. (flQ|) and (fTTjl . respectively. 




>4- 
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FIG. 9: Momentum distribution as a function of q in units of 
the lattice spacing. Note that at smaller U values there are 
three prominent interference peaks at q — and q — ±2tv. 
As U increases, the peaks get smaller, indicating a transition 
from SF to MI. 




FIG. 10: Condensate fraction (left) and visibility (right) as a 
function of U for a trap of depth Vt = 0.004. The inset zooms 
in on the kinks observed in the visibility corresponding to the 
formation of Mott shoulders. 



For a system in uniform SF phase, the fringe visibility 
is 1 (39j . In the homogeneous case, when the system un- 
dergoes a quantum phase transition from SF to MI, the 
visibility falls monotonously [29|, |4l| . In the presence of a 
trap, however, the visibility as a function of U has a rich 
structure due to the formation of alternating SF and MI 
shells [Hj]. From Figs. (JTOJ) and ([ITj). we note the follow- 
ing for the trap case: (1) the visibility remains finite even 
at high U values (i.e., deep inside the MI lobes) compared 




10 15 20 
U/t 



FIG. 11: Condensate fraction (left) and visibility (right) as a 
function of U for a trap of depth Vt = 0.008. 




FIG. 12: Variation of fio, the local chemical potential at the 
center of the trap as a function of U. For a given value of U, 
the inhomogeneous system can be represented by a vertical 
line originating at /iq. 



to the homogeneous case, (2) kinks develop over a cer- 
tain range of U and (3) the visibility drops drastically for 
particular values of U and for further increase in U, the 
variation is slow. Similar behavior is also noted in the 
condensate fraction. 

For a given value of U and number of bosons N, the ho- 
mogeneous system represents a point (//, U) in the phase 
diagram given in the Fig. ([6]). However, for an inhomo- 
geneous system the chemical potential varies across the 
lattice and is represented by a vertical line (see Fig. (fl~2"]) ), 
originating at fi for a given U in the (fj,—U) plane, where 
/xo is the chemical potential at the center of the trap. The 
values of /j,q as a function of U is shown in Fig. (fT2"|) . The 
behavior of the condensate fraction and visibility (shown 
in Figs. ([TU]) and (fTTj) ) is then easily understood by trac- 
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ing the po trajectory as a function of U. For U < 6.0, a 
vertical line starting at representing possible values 
of the local chemical potential for a given U, does not 
intersect the MI lobe and no Mott plateau forms in the 
density profile. As U increases, the system begins to fa- 
vor the MI phase, and as a result the condensate fraction 
and the visibility decreases monotonically. However, for 
U > 6.0, the vertical line, intersects the Mott lobe, re- 
sulting in a well-developed Mott plateau in the density 
profile. All the kinks in the condensate fraction and the 
visibility are observed for U > 6.0, indicating the forma- 
tion and broadening of the Mott plateau in the system, 
as the bosons re-distribute themselves between the two 
phases across the lattice. Finally for larger values of U, 
the fj,Q trajectory enters the Mott lobe and the central SF 
region vanishes completely. The entire system is in the 
Mott phase except for the edges. As a result, the conden- 
sate fraction and the visibility drops drastically. Further 
increase in U results in smooth variations of both these 
quantities as the SF phases exist only at the edges and 
the number of bosons in the SF phases do not vary much 
as shown in Fig. ©. 

In the experiments, the chemical potential is usually 
changed by changing the number of bosons N. We plot, 
in Fig. [)13p. the variation of the condensate fraction and 
the visibility as a function of N, for Vt — 0.008 and 
U = 10. We see that when the MI plateau forms, both 
visibility and condensate fraction dip, the first of these 
corresponding to the formation of a p = 1 Mott plateau, 
occurs around N ~ 40. This can be observed in the den- 
sity profile (see Fig. [fTl))). The plateau in Fig. [[15]). in the 
condensate fraction, indicates the range of N for which 
the p = 1 MI phase is possible for the given U value. We 
note that beyond N ~ 56, a 1 < p < 2 SF forms, marked 
by an increase in the visibility and the condensate frac- 
tion. The second minimum occurs around N ~ 140 sig- 
nals the formation of the second Mott plateau (p = 2), 
as shown in Fig. (fl"4|) . Further increase in N, beyond 
N ~ 150, results in another SF phase, with the on-site 
density ranging between 2 < p < 3. Therefore, for a 
fixed value of U, the minima in the condensate fraction 
and the visibility as a function of N are good indicators 
for the formation of Mott plateaus. 




FIG. 13: Condensate fraction (left) and visibility (right) as a 
function of N the total number of bosons in a trap of depth 
Vt — 0.008 and U — 10. Both the condensate fraction and 
visibility show local minima for particular ranges of N, indi- 
cating the formation of MI phases. The subsequent increase 
in these quantities indicate the formation of SF phases. 
3 



N = 30 
N = 44 



N = 60 

N = 100 

N = 140 

N= 150 




FIG. 14: Density profile for different total number of bosons 
N. Note that the minima in visibility and condensate fraction 
in Fig. (|13|l correspond to formation of MI plateaus. 



V. CONCLUSION 

This paper demonstrates a way of extracting the phase 
diagram for the Bose-Hubbard model in the presence of 
an external trapping potential, from the number density 
profile. Signatures of these phases in experimentally ob- 
served quantities such as visibility, condensate fraction 
and N(p) have been documented. We have also obtained 
the density-profile as a function of the chemical potential 
for the homogeneous case, using FS-DMRG for the first 
time, to the best of our knowledge. 

Future directions are immense. As a first step, the 
extended Bose-Hubbard model can be investigated and 



the phase diagram, together with the experimental sig- 
natures for the various phases can be extracted in an 
analogous way, which is in progress. The influence of 
a three-body term for the Bose-Hubbard and the ex- 
tended Bose-Hubbard can be investigated, which is also 
in progress. 

The presence of a trap makes the system interesting 
due to the simultaneous existence of different phases 
(such as the SF and the MI phase in this work) and 
gets the theoretical predictions closer to experiments. 
It would also be interesting to study the Bose-Hubbard 
model for two boson species and Bose-Fermi mixtures in 
the presence of a trap. This paper serves more as bench- 
mark to extract the phase diagram in a straight-forward 
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